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Abstract 

This paper concerns spectral stability and time evolution of 7V-solitons in the KdV hierarchy 
with mixed commuting time flows. Spectral stability problem is analyzed by using a pair of 
self-adjoint operators with finite numbers of negative eigenvalues. We show that the absence of 
unstable eigenvalues in the stability problem is related to the absence of negative eigenvalues 
of these operators in the constrained function spaces. Time evolution of ./V-solitons is uniquely 
characterized from the inverse scattering transform technique. 



1 Introduction 

We address the stability problem for solitary waves in the KdV-type evolution equation: 

du d 5H ^ ^ 

dt dx 5u ' 

where u G K, x G M, t G M+, and H(u) is the Hamiltonian. Besides the canonical Korteweg-de 
Vries (KdV) equation |N85j , the time-evolution problem Ql.l|) includes higher-order KdV equations 
of the integrable KdV hierarchy (e.g. see |KJN781 HTT78 ). The integrable KdV hierarchy can be 
derived with the asymptotic multi-scale expansion technique for modelling of solitary waves in 
physical non-integrable problems, e.g. in ferromag nets [TIM lL03j (see also |HK02| for the normal 
form of perturbed KdV equation). 

Our work originates from analysis of spectral stability of solitary waves in the KdV-type evolution 
equations BSS87[ rW87|ISS9C)|lPW 92_. Let us assume that the evolution problem 1)1. lj) has a solitary 
wave solution u = uo(x,t) that decays exponentially in space x G K and changes in time t G 1R+ 
according to symmetries of The spectral stability problem for KdV solitary waves takes the 

general form: 

d x Cv = Xv, (1.2) 

where v £ C, A G C, and the self-adjoint linearized operator C is computed at the solution uo(x,t) 
after separation of variables (x,t), i.e. C = D 2 H{uq). The solitary wave solution uo(x,t) is 
spectrally unstable if there exist an eigenvalue A G C with Re(A) > and v G L 2 (IR). It is weakly 
spectrally stable if no eigenvalues A G C with Re(A) > and v G L 2 (W) exist. 
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We shall characterize unstable eigenvalues of the stability problem (|1.2|) from the study of the 
self-adjoint operator C. The operator C defines the energy quadratic form: 

h(v) = (v,Cv) := / v{x){Cv){x)dx. (1.3) 
Jr 

The relation between eigenvalues of the stability problem and those of the linearized energy was 
recently studied in the context of spectral stability of solitary waves in the nonlinear Schrodinger 
equations C PV051 IP05| . The number of unstable eigenvalues in the spectral stability problem 
was found to be related to the number of negative eigenvalues of the energy quadratic form in a 
constrained function space. The same relation was derived independently in KKSU 11 by extending 
earlier results |G90| to an abstract Hamiltonian dynamical system. 

The methods of KKS04, IP05| are not applicable to the KdV-type evolution equations, since 
the symplectic operator d x is not invertible. Nevertheless, we will show with explicit analysis of 
quadratic forms that the spectral stability problem (|1.2[) can be embedded into a larger problem, 
which has the same structure as that considered in KKSfiH IP05| . Since the operator C maps L 2 (R) 
to L 2 (R), the eigenfunction v (x) G L 2 (R) for A ^ must satisfy the constraint: 

(l,v) = / v(x)dx = 0. (1.4) 
Jr 

If v G L 2 (R) decays exponentially as \x\ — ► oo and satisfies the constraint ()1.4|) . the eigenfunction 
of the spectral problem (|1.2|) can be represented as v = —w'(x), where w G L 2 (R), i.e. v G im(<9 x .). 
We shall hence replace the stability problem (|1.2|) by the coupled system: 

Cv = -Xw, Mw = Xv, (1.5) 

where M = —d x Cd x . Eliminating w from the system Q1.5JI . we reduce the coupled problem ()1.5|) 
to the scalar problem: 

d x Cd x Cv = X 2 v, (1.6) 

which has both eigenvalues A and —A of the original problem (jl.2|) . Moreover, if A is a simple 
eigenvalue of d x C, then the system (|1.5j) has the solution v = —w'(x), while if A is a simple 
eigenvalue of —d x C, then the system ljl.5|) has the solution v = w'(x), 

The coupled system (jl . 5|) is defined in a constrained subspace of L 2 (M, C 2 ). Assuming that 
the linear operator C has an isolated kernel and a positive continuous spectrum and that the 
eigenfunctions of ker(£) satisfy the constraint (|1.4|) . we introduce the constrained subspaces: 

X C (R) = {v G L 2 (R) : («, ker^)) = 0} , (1.7) 
X' C (R) = {w G L 2 (R) : (w,ker(£)) = 0} . (1.8) 

If (v, w) G L 2 (R, C 2 ) is the eigenvector of the system (fT3|) for A / 0, then v G X C (R) and w G X' C (R). 

Within this general formalism, we shall study A^-solitons of the KdV hierarchy. Time evolution 
of the iV-solitons is defined by the mixed commuting flows of the higher-order KdV equations. 
Using earlier results of MS93 , we shall prove that the operators £ and Ai have no negative 
eigenvalues in the constrained spaces X C (IR) and A^(M), respectively. This fact has a consequence 
that the stability problem (|1.2|) and the coupled system (|1.5|) have no unstable eigenvalues with 
Re(A) > 0. We also use the inverse scattering transform technique AK82] to characterize the 
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time evolution of iV-solitons in the mixed commuting flows of the KdV hierarchy. Uniqueness of 
./V-solitons in a constrained variational problem is also proved. Both main results (spectral stability 
and time evolution) are not covered by the previous publications on iV-solitons in the KdV hierarchy 
|MS93MMMT02| . 

The paper is organized as follows. In Section 2, we briefly summarize the background of the 
integrable KdV hierarchy, and give an explicit definition of the spaces X C (K.) and -X"£(IR) for the 
iV-solitons. Spectral stability of iV-solitons in mixed commuting time flows is studied in Section 3. 
Time evolution of iV-solitons is characterized in Section 4. Section 5 concludes the paper. 



2 Review of iV-solitons in the KdV hierarchy 



The integrable KdV hierarchy (see review in |N85| ) is defined by the set of KdV- type evolution 
equations: 

du d 5H n 

df n = Tx^ neN+ ' (2 ' 1} 

where H n (u) are the Hamiltonian in energy space fP -1 (R). The Hamiltonians are constructed 
recursively as 

5H n+ i 5H n 

J — ~ = K.— — , n G N_|_, (2.2) 

ou ou 

where the linear operators J and JC take the form: 

d d 3 ( d d\ 

ox ox- 3 \ox Ox J 

and the lowest three Hamiltonians are: 

H x = \ [ u 2 dx, (2.4) 

z JR 

H 2 = \!{ul- 2u 3 ) dx, (2.5) 

H 3 = - (u 2 xx - Wuul + bu 4 ) dx. (2.6) 

The Hamiltonians H\, H2, and H% generate the first three members of the KdV hierarchy, which 
are given by the transport equation = u x , the KdV equation, 

u* 2 = -u xxx - 6uu x , (2.7) 

and the integrable fifth-order KdV equation, 

ut s = u xxxxx + Wuu xxx + 20u x u xx + 30u 2 u x . (2-8) 

All members of the KdV hierachy (|2.1|l - ()2.6jl have families of (2iV)-parameter solutions called 
iV-solitons which are expressed by the r- function |N85j . 

d 2 

u(x;t 1 ,t 2 ,...) = U N (x - 6»i, ...,x - 9 N ) = 2^ logr(x - U ...,x - 9 N ), (2.9) 
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where 



n (ti, t 2 , ...) = S n + ^(-l) fc c£ x 4 = 5 n - ti + c n t 2 - 4*3 + 



(2.10) 



fc=i 



and n = 1,2, ...,N. Parameters <5 n are arbitrary, while parameters c n satisfy the constraints c n > 



for all n = 1,2,..., N. 



We give explicit formulae for the first two solitons ()2.9|) ~ (|2.1(J[) . The 1-soliton is given by: 



d 2 

Ux = 2— log 

ox z 



1 _|_ e Vci( x - d i) 



ci > 0. 



The function U\{x) is a decaying solution of the second-order ODE: 



H' 2 (u) + v 1 H , 1 (u) 
where v\ = c\. The 2-solitons are given by: 



-u xx - 3u + c\u = 0, 



d 2 

U 2 = 2— log 

ox z 



1 _|_ gv 7 ^"^- 61 !) _)_ e V^(x-9t) _|_ 



,V^T(a!-0i)4-v/ra(a!-^2) 



(2.11) 
(2.12) 

(2.13) 



(2.14) 



where c\, c 2 > and ci ^ c%. The function ^(x) is a decaying solution of the fourth-order ODE: 

H'z(u) + v 2 H' 2 {u) + Vl H[(u) = 
u xxxx + 10uu xx + 5ul + 10u 3 - (CX + c 2 )(u xx + 3u 2 ) + cic 2 u = 0, 

where v\ = c\c 2 and v 2 = c\ + c 2 . 

In general, the functions Un(x) are critical points of the Lyapunov functional in H N ( 

N 

A N (u) = H N+1 (u) + v n H n (u), 



(2.15) 



such that 



n=l 



A? 



A' N {U N ) = H' N+1 (U N ) + £ u n H' n {U N ) = 0, 



(2.16) 



71=1 



where the Lagrange multipliers v\,...,vj^ are elementary symmetric functions of ci,...,cjv due to nor- 
malization of H n (see |MS93 ). Lyapunov stability of iV-solitons as critical points of the constrained 
Hamiltonian Hn+i(u) was proved in |MS93j . Specifically, 7V-solitons Un(x) are minimal points of 
Hn+i(u) subject to N constraints on the lower-order Hamiltonians: 



H n {u) = constant, n = 1, N, 
so that the second variation of A^(u) is positive definite, 



(v,C N v) 



hm ^ > 0, v G X c 



e^0 



nx' 



(2.17) 



(2.18) 



Here Cn is the self-adjoint linearized operator of (2A r )-order with finite number of negative eigen- 
values, finite-dimensional kernel and positive continuous spectrum, bounded away from zero (see 
MS93 ), -Xc(R) is the closed orthogonal compliment of the kernel of Cn8 x : 



X c 



v G L 2 (R) : v, 



SU N J 



n=l,...,N\, 



(2.19) 
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and X' C (M) is the closed orthogonal compliment of the kernel of Cm- 

X'M) = [™ G L\R) : =0, n = l,...,iv}. (2.20) 



The embedding u € X C (R) follows from the constraints (|2.17j> . while the embedding v £ -X^ 
is set artificially, in order to remove the zero eigenvalues of Cm and to ensure positivity of the 
energy quadratic form (|2.18|) . Due to positivity ([2.18ft . the functional A(it) is convex at the point 
u = Um{x), such that Lyapunov stability of iV-solitons in energy space H N (M) holds _MS93 . 

The linearized operators C\ and C 2 for 1-soliton and 2-solitons are given explicitly as: 

£i = -a£-6C/i(x) + ci, (2.21) 

and 

C 2 = 8i + \W 2 {x)d 2 x + lW 2 {x)d x + 10t#(z) + 30CT|(x) - (ci + c 2 ) (fig + 6U 2 {x)) + cic 2 . (2.22) 

The first two eigenfunctions of the kernel of Cm are: 

Vl = U' N (x), v 2 = -U%(x) - 6U N (x)U' N (x). (2.23) 

The first two eigenfunctions of the kernel of Cm9 x are: 

Wl = U N {x), w 2 = -U N (x) - 3U N (x). (2.24) 

These explicit expressions illustrate the general construction of the linearized operator Cm and the 
kernels of Cm and CmO x . 



3 Spectral stability of iV-solitons 

Time evolution of iV-solitons ([2.9|) ~ ([2.1U|) is defined in the mixed commuting flows of the KdV 
hierarchy, which are generated by the Lyapunov functional Q2.15|) : 

du du A du d 5Am , , 

dt dtM+i " dt n dx 5u 

n=l 

The family of iV-solitons u = Um(x; ci, cm', 5%, 6m) is a time-independent space-decaying solu- 
tion of the KdV-type evolution equation (J3.1|) . where ui,...,pm are elementary symmetric functions 
of ci,...,cjv- Linearization of the evolution equation 1)3. 1|) as u(x;t) = Un(x) + Vn(x; t) and separa- 
tion of variables as Vjv(x;i) = v(x)e^ t results in the spectral stability problem 

8 x Cmv = Xv (3.2) 

where Cm is the same as in the energy quadratic form (|2.18|) . We assume that the family of N- 
solitons is non-degenerate MS93 and characterize the kernel of Cm in terms of the symmetries of 
the KdV hierarchy |KT78j . 

Assumption 3.1 N-solitons Um(x;ci, ...,cm',5i, ...,5n) have distinct positive parameters c\,..., cm ■ 
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Lemma 3.2 The kernel of Cjy has a basis of N linearly independent eigenf unctions {v n (x)}^ =1 in 
L 2 (R), where 

v n (x) = -f^, n = l,...,N. (3.3) 

OX OU TV 

Proof. Derivatives of Un{x; c\, c/y; S\, $n) with respect to arbitrary parameters 5\,...,5n form 
a basis of the kernel of £at, since they are linearly independent and Cn is a self-adjoint operator 
of (2iV)-order. It follows from (j2~TUj) that 

f^B-ir'c^B-ir'cr^, B =i,...,*. (3 .4, 

fe=i fc=i 
Since the Vandermonde determinant of ci,...,c/v is non-singular under Assumption l3.il the basis of 
{fsf } x is ec l uivalent to tlie basis {^f } • Jt follows from the KdV hierarchy ({2~T|) that 

See also |KT78j and |MS9 3, Lemma 3.4] for alternative proofs. □ 



Lemma 3.3 The generalized kernel of d x C]\f has a basis of N linearly independent eigenf unctions 
{v { n\x)}% =1 inL 2 (M), where 

v^(x) = -^, n = l,...,N. (3.6) 

Proof. The generalized kernel of d x C n is generated by solutions of the inhomogeneous problem: 

8 x C N v^ =v n (x). (3.7) 

Integrating Q3.7JI in x for Vn^ £ L 2 (M), we find that eigenfunctions Vn\x) satisfy the inhomogeneous 
equations, 

C N v^ = ^, n = l,...,N. (3.8) 
dU N 

It follows from the derivative of the variation equations (|2.16j) in u n that 

8Un SH n 

C N — — = -— — , n = l,...,N. 3.9) 
du n 5U N 

As a result, relations ()3.6|) hold. □ 



Remark 3.4 The symmetries v n (x) and Vn (x) are shown to be expressed by the squared eigen- 
functions of the inverse spectral method IKN78\ \KT78\[ . The squared eigenfunctions form a basis 
for a class of functions u(x) satisfying the bounds: 

/oo 
(1 + x 2 )\u(x)\ dx < 00. (3.10) 



6 



Define the solution surface as A s = Ajy(U]y) = A s (ui, v^). The Hessian matrix H of principal 
curvatures of the solution surface A s (^i, vn) has the elements: 

d 2 A s dH n (U N ) dH m {U N ) 

Ti-n,m = -^ — 7\ — = — ^ = a ' n,m = 1, N, (3.11) 

where 

dH n {U N ) (8U N 5H n \ (8U N du N \ 

— a = I "a — >I7T~ l = ~ « — > Cn ^ — ' n,m = 1, N. (3.12) 

dv m \ou m dU N ) \ov m du n J 

Let z{TL) be the number of zero eigenvalues of 7i. 

Lemma 3.5 The second generalized kernel of 8 x Cn is empty in L 2 (R) if z(TL) = 0. 

Proof. The second generalized kernel is generated by solutions of the inhomogeneous problem: 

N 

8 x C N v^ = Y j anv£ ) {x), (3.13) 

n=l 

where {a n }^ =1 is a set of coefficients. Computing the inner products of the inhomogeneous problem 
(|3.13j) with the set of eigenfunctions {|^ IL }^ r =1 , we find that the solutions v^ 2 ' G L 2 (R) may exist 
only if the Hessian matrix TC in (|3.11|) - (|3^12|) has a zero eigenvalue. They do not exist if z(TC) = 0. 
□ 



Lemma 3.6 The kernel of A4n = —d x CNd x has a basis of N linearly independent eigenfunctions 
{w n (x)}n=i in L 2 (M.), where 

5H 

w n (x) = —^, n = l,...,N, (3.14) 
dU N 

Proof. Integrating d x £Nd x w = in x, we have Cnw'{x) = for w'(x) G L 2 (M). Therefore, 
w' G ker(£7v)j such that the eigenfunctions (|3.14f) follow from integration of eigenfunctions (|3.3[) . 
□ 



We now work with the coupled problem 

Cnv = —Xw, M.nw = Xv, (3.15) 

which is equivalent to the spectral problem 1)3. 2 Jl when v = —w'(x), subject to the constraint 1)1.4)1 . 
By Lemmas 13. 21 and 13 . 6) definitions 1)1.7)1 - 1)1.8)) of the constrained spaces X C (M) and X' C (M) coincide 
with the definitions (|2~T§|) - l(2~2lH) . It is therefore clear that v G X C (R) and w G X' C (R) for any 
solution (v,w) G L 2 (R, C 2 ) of the coupled problem (|3~T5|) with A / 0. 

Analysis of unstable eigenvalues in the coupled problem (|3.15|) is based on the fact that both 
operators Cn and Mn have finitely many negative eigenvalues in L 2 (R). The continuous spectrum 
of Cn is positive and bounded away from zero by Co = n^=i °n > 0, due to the exponential decay 
of potential functions and the factorization relation MS93 : 

U N {x) = 0: C N = (-d 2 + a)...(-dl + c N ). (3.16) 
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The continuous spectrum of M.n is non-negative, due to the explicit relation: M.n = —d x £-Nd x - 
The kernel of Cn belongs to X C (R), due to commutativity of the Hamiltonians H n (u) of the KdV 
hierarchy $n $ -$Z ty : 

(5H n d 5H m \ 

Equivalently, the kernel of Mn belongs to X' C (M). These facts enable us to consider the number of 
negative eigenvalues of Cn and M.^ in X C (NL) and X' C (M.), respectively, and to relate the absence 
of negative eigenvalues to spectral stability of iV-solitons in the spectral problem (|3.2|) . 

Proposition 3.7 Let u{Cn) be the number of negative eigenvalues of Cn in L 2 (R) andp{H) be the 
number of positive eigenvalues ofH. The numbers #<o(£jv) and #<o{M.n) of negative eigenvalues 
of Cn in X C (1R) and of M.n in X' C (M.), respectively, are given by: 

#<o(£n) = #<o(M N ) = n{C N ) -p{H). 



Proof. It follows from jSEESOl Theorem 3.1] and |MS93I Lemma 2.1] that #<o{£ N ) = n(£ N )-p(H). 
Equivalently, the same result can be proved by minimization of quadratic forms with Lagrange 
multipliers, see PQ5, ICPV05) . In order to prove the relation #<o(-Mat) = #<o(£iv)> we integrate 
the quadratic form (u>, M.j^w), w G X' C (M) by parts: 

(w, M-nw) = — (w,d x £Ndxw) = (w',£nw') ■ (3.18) 

By integrating the constraints in the inner products (|2,2U|) by parts, we confirm that if w G X' C (M), 
then w' G X C (R), such that the quadratic form (w',Cnw') is defined in w' G X C (R). Finally, it 
follows from (|3.3() and (|3.14|) after integration by parts that 

{vn^M^Vm) = (w n ,C~^~w m ) , n,m = l,...,N, (3.19) 

where the functions M.^ v m are bounded due to the fact that v m {x) satisfy the constraint (|1.4j) for 
any m. Applying Lemma 3.4 from jCPV05| . we have the relation #<o(A^at) = ti(Cn) — p(7i) = 
#<o(Cn). □ 



Proposition 3.8 Let p(Tl) = n(£jv). The coupled problem U.5\) has no eigenvalues A G C with 
Re (A) > and (v,w) G L 2 (M,C 2 ). 

Proof. When p{TL) = n(£jy), we have by Proposition 13.71 

Vu G I c (l) : (v,C N v)>0, (3.20) 
VweX' c (R) : (w,M N w) > 0, (3.21) 

where the zero values are achieved if and only if v G ker(£jv) i n -Xc(R) an d w G ker(.M n) i n -^c(^)- 
We assume that there exists an eigenvalue Ao G C with Re(Ao) > and (vo,wo) G L 2 (R, C 2 ) and 
show the contradiction. If the eigenvalue Ao would exist, then the coupled system (|3.15|) would 
result in the identity: 

Xq(w q , M n w ) + \ (vq,Cnvo) = 0, (3.22) 
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where the real- valued inner products are used for clarity of notations. Since Re(Ao) > and 
quadratic forms associated with Cn and M.n are real- valued, we have 

Oo, M N w ) = -(v ,£nv ), 

which is the contradiction, since both quadratic forms are non-negative for vq E X c (M) and wo 6 
X' C (M.) and non-zero for vq £ ker(£/v) and wq ^ ker(A^Ar) (when Aq 7^ 0). □ 



Spectral Stability Theorem: Let N-solitons Un(x) of the KdV hierarchy \2.1\) - \2.b}) satisfy 
Assumvtion \3.1\ Then, N-solitons are weakly spectrally stable such that the spectral problem h3. 2]) 
has no eigenvalues A G C with Re(A) > and v G L 2 (R, C). 



Proof. We show that conditions of Proposition 13.81 are satisfied for A^-solitons of the KdV hierarchy 
(12TT |H |21) |) . It was proved in |MS931 Lemma 3.5] and in )MS931 Lemma 3.6] that 



n(C N ) 



N + 1 



p(H) 



N + 1 



(3.23) 



where [z] is the integer part of z. Therefore, the result holds by Proposition 



□ 



Remark 3.9 It was also proved in 'M S93\ Lemma 3.7] that z(H) = 0. Therefore, the second 
generalized kernel of 8 x Cn is empty by Lemma \ °3.5l such that the zero eigenvalue A = in the 
spectral problem \1.2\) is controlled by the symmetries of the KdV hierarchy. 



4 Time evolution of 7V-solitons in the KdV hierarchy 

According to l|2.18|) . the Lyapunov functional Ajv(ii) defined in (|2.15|) is convex at the point u = 
Un(x), where Un(x) is iV-solitons of the KdV hierarchy (|2.1)) - (|2.6|) . We give a direct proof of 
convexity based on the inverse spectral method |N85j . Moreover, we prove uniqueness of iV-solitons 
as minimizers of the constrained variational problem 1)2. 16(1 - 1)2. 17j) . Using the asymptotic results 
from AK82 , we formulate the Time Evolution Theorem for iV-solitons in the energy space H (R) , 
which follows from Lyapunov stability of A^-solitons. 

These results improve the Lyapunov Stability Theorem of MS93 , where the time evolution of 
parameters n (t) in the solution Un(x) remains undefined and uniqueness of minimizers in the 
constrained variational problem ()2. 16(1 - 1)2. 17(1 remains open. The unique characterization of the 
time evolution of parameters of N solitons is equivalent to the asymptotic stability of solitons in 
the mixed commuting time flows ()3.1(l . We note that asymptotic stability of 1-soliton in the KdV 
equation was proved in )PW94j with exponentially weighted spaces (see also |B72|IB75] for pioneer 
papers), while asymptotic stability of A-solitons in the energy space H l {K) was recently proved in 
MMT02) Corollary 1] with different analysis of energy functionals. 

Let us recall some necessary formulae obtained in the inverse spectral method N85 . The method 
is based on the isomorphism between the potential u(x, t) of the Schrodinger equation and the 
scattering data S(t) (see, e.g., |N85j ) . Consider the scattering problem, 

(dl + u(x,t) + k 2 ) V = 0, (4.1) 
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with the boundary conditions, 

V>(x,t;A;) — ► | a (Q e -ikx + b ( t .fy e ikx as x ^oo ■ ( 42 ) 

Assuming that the potential u(x, t) is exponentially decaying in x, it is proved that the function 
a(k) is analytic on the upper half plane of k £ C, and the bound states are defined by the zeros of 
a(k) at k = tKj, Kj > 0, such that 

a{%Kj) = j = l,...,N. (4.3) 

The scattering data is then defined by 

S(t) = [{K 3 ,C,(t)}f =1 , {r(t;k)} km ] , (4.4) 

where r(t; k) = b(t; k)/a(k) is the reflection coefficient, and Cj(t) is the normalization constant (see 
N85 for details). The parameter Cj in the iV-solitons (|2.9|) - (|2.1U[) is given by Cj = Ak 2 - for all j. 
We assume that the parameters {%}|Li are ordered as 

Ki > K 2 > ■ ■ ■ > KN > 0. (4.5) 

By inverse scattering, the potential u(x, t) can be expressed as 

N 2i r°° 

u(x,t) = A s ^K j C j {t)tp 2 {x,t;iK j ) H / kr(t;k)i(j 2 (x,t;k) dk , (4.6) 

i=i " " " ' x 

which shows that u(x,t) consists of A^-solitons and the radiation. With the formula (|4.6j) . the 
Hamiltonian H n can be expressed in terms of the scattering data N85 , 



H n {u) = (-Ijn+lj^ 1 f^VlE K f +1 " (-!) n ^j 



(4.7) 



where R n is the radiation part, given by 



1 



oo 



R n = — - / A; 2n ln(l - \r(t; k)\ 2 ) dk > 0, n=l,2,... (4.8) 
27r Jo 

We show that the alternate sign (—1)" in front of R n in the representation H4.7|) plays a crucial rule 
for the convexity. 

Proposition 4.1 The N-soliton Un(x) is uniquely determined in the variational problem \2. 1 6}) 
by the constraints \2.17\l except the phase parameters S±, . . . , 5jy. The Lyapunov functional An(u) 
is convex at u = Un(x). 

Proof. Let us consider the variation of H n (v) at u = Un(x), which we denote as Ai? n (£/jv)- The 
constraints ()2.17|) imply that AH n (UN) = 0, such that 

N 

^2K 2n A Kj - {-l) n AR n = 0, n = l,...,N, (4.9) 
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where Akj is the variation of Kj, and AR n is the variation of R n . Since R n = at u = Ujy(x), we 
understand from the explicit formula (|4.8j) that all AR n > 0. It is clear from the system (|4.9|) that 
the variations of Kj are balanced with the variations of R n , so that the Hamiltonian H n remains 
constant. The variations AUn belongs then to the constrained space X C (R), defined in (|2.19j) . We 
also note that the number of solitons may change under the variation of H n (u), but this does not 
appear in the constraints Q4.9JI . since Kj = for j > N. We need to prove the convexity of the 
functional Ajv(u) at u = U]\r(x), which implies that 

AH N+l (U N ) > 0. (4.10) 

^From the system (|4.9j) . the variation AHn + i(Un) can be expressed as 

N 

AH N+1 (U N ) = (-l) N 2 2N+3 [S2Kf^AK, + (-l) N AR N+1 



^ ( DiARi + (-l) N ARf 



N 



(4.11) 



vZ=l 



where D is the determinant of the N x N matrix Kn '■= {{^) '■ 1 < hj < ^} an d Di is the 
determinant of the matrix Kn after replacing the l-th. row by {{k 2N+2 ) : 1 < i < N}. It is clear 
that D is the Vandermonde determinant, computed as 

D= H (4.12) 

l<i<j<N 

so that -D > from the ordering (|4.5[) . We further notice that Di can be written as 

A = (-l)^^^ 1 , (4.13) 

where £>f +1 is the determinant of the (iV + 1) X (iV + 1) matrix Kjv+i = {{k 2j ) : 1 < i,j < N + l} 
after removing the Z-th row and the (N + l)-th column. Thus we have 

02JV+2 / N \ 
AH N+1 (U N ) = — — \Y2D? +1 AR l + Ai?Ar +1 J . (4.14) 

Since all Ai? n > 0, we only need to show that -D/ V+1 > for the convexity. This follows from the 
expression, 

D? +1 = a N a N ^ l+1 D, (4.15) 
where Ok is the symmetric polynomial of (k 2 , . . . , k 2 n ) given by 

a k = y, 4 •••4 >°- ( 4 - 16 ) 

l<h< — <i k <N 

As a result, D^ +1 > and the convexity (|4,1U|) holds if there exists at least one AR n ^ 0. (Note 
that if ARn ^ for some n, then it is not zero for any n.) Otherwise, i.e. when all ARn = 0, 
all Ak h = from the system (|4.9|) . and the variations of Un(x) simply translate the solution 
Un(x) along its phase parameters 5i,...,5n- When the variation AUjy G X C (M) belongs also to the 
constrained space X' C (M), defined in 1)2.20(1 . the latter case is excluded. Thus the minimal point 
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of Hjy + i is given by the iV-soliton u = Un(x), that is, all R n = 0. Uniqueness of the minimizer 
Un(x) follows from the geometry. Each constraint H n (u) = constant gives a hypersurface in 
K2, «at) £ • If all those hypersurfaces intersect transversally with nonempty intersection, 
we have N\ points of the intersection. Then the ordering k\ > ■ ■ ■ > kn > chooses a unique point 
of the intersection. □ 



Time Evolution Theorem: For all e > 0, there exists 5 > such that if \\u(x, 0)— Un(x)\\ h nw\ < 
5, then there exist dp{t) = Cjt + Sj, where t = ti, such that 

\\u(x,t)-U N (x-0f(t),...,x-9x(t))\\ L 2(\ x \ >CNm <e, as t -> ±oo. (4.17) 

Proof. The proof relies on Proposition 14. II and formalism above. The time evolution of iV-solitons 
has been obtained in |AK82j by the inverse spectral method, see Eqs. (2.15a) and (2.15b) in 
AK82 a , where n^t should be read as An 2 N t = cjyt. We only note that the soliton speeds Cj = 4k^ 
and soliton phases Sj must be obtained by the scattering problem (|4.1|) for the perturbed potential 
u = U N + AU N . □ 



5 Summary 

We have studied the spectral stability and time evolution of iV-solitons in the mixed commuting 
flows of the integrable KdV hierarchy. We have proved that iV-solitons of the KdV hierarchy 
are spectrally stable with respect to the time evolution. The analysis extend the recent results of 
CPV05, KKS04, P05 to the KdV- type evolution equation Ql.l|) . In particular, the proof of spectral 
stability is related to the absence of negative eigenvalues of linearized energy in constrained function 
spaces. We have also characterized time evolution of iV-solitons and proved that the iV-solitons 
are uniquely determined in a constrained variational problem. Further development of the spectral 
stability analysis to non-integrable KdV-type evolution equations, e.g. to the fifth-order KdV 
model, will be reported elsewhere. 

Acknowledgements: Y.K. is partially supported by the NSF grant DMS-0404931. D.P. is 
partially supported by the NSERC Discovery grant. 
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